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, We consider Dirac fermions moving in a plane with a static homogeneous magnetic field or- 

thogonal to the plane. We calculate the effective action at finite temperature and density. The 
magnetization is derived and it is shown that the fermion gas exhibits de Haas-van Alphen oscil- 
lations at small temperatures and weak magnetic fields. We also comment upon earlier work. 
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1 Introduction 

a; 

Many of the phenomena that have been discovered in condensed matter physics over 
the last few decades are to a very good approximation two dimensional. The most 
important of these are the (Fractional) Quantum Hall effect and high T c supercon- 
ductivity pj. 
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Quantum field theories in lower dimensions have therefore become of increasing 
interest in recent years. Both systems mentioned above have been modelled by 
anyons, which are particles or excitations that obey fractional statistics. Anyons 
can be described in terms of Chern-Simons field theories [1-3]. The literature on 
2+ld field theory in general and Chern-Simons field theory in particular is now vast 
(see e.g refs. [1-5]), and we shall only comment upon a few selected papers, which 
are relevant for the present work. 

Ten years ago Redlich |J considered fermions in a plane moving in a constant 
electromagnetic field. Using Schwinger's proper time method to obtain the ef- 
fective action for the gauge field, he demonstrated that a Chern-Simons term was 
induced by radiative corrections. The Chern-Simons term is parity breaking and 
is gauge-invariant modulo surface terms. More recently, Lykken et al. have ex- 
amined two dimensional Dirac fermions coupled to a gauge field whose dynamics 



1 



is governed by a (topological) Chern- Simons term. It was found that this planar 
Fermi gas becomes superconducting (below a critical temperature) when the induced 
Chern-Simons term exactly cancels the one appearing in the classical action. 

External electromagnetic fields may give rise to induced charges in the Dirac 
vacuum if the energy spectrum is asymmetric with respect to some arbitrarily chosen 
zero point. The vacuum charge comes about since the number of particles gets 
reduced (or increased) relative to the free case. Furthermore, induced currents may 
appear and are attributed to the drift of the induced charges. This only happens if 
the external field does not respect the translational symmetry of the system. These 
interesting phenomena have been examined in detail by Flekk0y and Leinaas || in 
connection with magnetic vortices and their relevance to the Hall effect has been 
studied by Fumita and Shizuya 0. 

In the present paper we re-examine the system considered by Redlich M. We 
shall restrict ourselves to the case of a constant magnetic field, but we extend the 
analysis by including thermal effects and we shall mainly focus on the magnetization 
of the system. 

In the first section we find the solutions of the Dirac equation for particles in 
a constant magnetic field. These solutions are then used to calculate the thermal 
fermion propagator, from which we derive the effective action for the gauge field. In 
the fourth section we compute the magnetization. It is shown that the de Haas-van 
Alphen oscillations [10,11] occur at low temperatures and weak magnetic fields, and 
that is the main result of this paper. We consider various limits of the magnetization 
that may be understood from a physical point of view and we explicitly demonstrate 
that we get the correct high temperature limit (i.e. free fermions). 

Finally, we summarize and draw some conclusions in the last section. 



2 Dirac Equation 

In this section we shall briefly discuss the solutions of the Dirac equation in 2+ld 
with a constant magnetic field along the z-axis. The Dirac equation reads 

(i<y»d li -erfA lt -rn)il> = Q, (1) 

where the gamma matrices satisfy the Clifford algebra 

{Y,Y} = 2g^. (2) 

In 2+1 dimensions the fundamental representation of the Clifford algebra is given by 
2x2 matrices and these can be constructed from the Pauli matrices. Furthermore, in 



2 



2+ld there are two inequivivalent choices of the gamma matrices, which corresponds 
to 7 M — > — 7 M . From eq. (ED) we see that this extra degree of freedom may be 
absorbed in the sign of m. These choices correspond to "spin up" and "spin down", 
respectively. 

We would also like to point out that the angular momentum operator a is a 
pseudo vector in 2+ld, implying that the Dirac equation written in terms of these 
matrices does not respect parity. This is no longer the case if the Dirac equation 
is expressed in terms of 4 x 4 matrices. The latter representation is reducible and 
reduces to the two inequivalent fundamental representations mentioned above. Let 
us explicitly verify this. A natural choice is the three block diagonal 4x4 matrices 




where 7 , 7 1 and 7 2 themselves are 2x2 matrices satisfying the algebra. This 
representation trivially reduces to the sum of the two inequivalent choices of the 
2x2 gamma matrices. Furthermore, the Dirac equation written in terms of the 
above matrices respects parity (a parity transformation simply interchanges the 
coupled equations and thereby the components of the spinor). 

In the following we make the choice 7 = 03, 7 1 = —i<72, 7 2 = — i&i and we have 
chosen the asymmetric gauge = (0, 0, —Bx). The metric is diag (—1, 1, 1). 

The Dirac equation then takes the form 



rn — ieBx + 

—i-^ — ieBx + — i-r^ — m 



^(x,t) = 0. (4) 



Here ip(x.,t) is a two component spinor. We now assume that the wavefunctions 
may be written as 

Vfc(x, t) = exp(-zEt + iky) ( ^ j , (5) 

where k denotes all quantum numbers necessary in order to completely characterize 
the solutions . Inserting this into eq. (SI) one obtains 



E — m — £ + 

—E — m 




(6) 



where 

t ± = -id x Ti{k-eBx). (7) 



3 



The equation for f K (x) is readily found from eq. (@): 

- m 2 - f K (x) = 0. 



The eigenfunctions of provided that e£> > 0, are JT2 

fc 



i»,jfc = ( — ) 4 exp 

7T 



\FleB{x- — ) 
ei3 



(9) 



where H n (x) is the n'th Hermite polynomial. Furthermore, I n ,k{x) is normalized to 
unity and satisfies 



0*0 = -W2eBnI n ^ k (x), 
£+I n ,k(x) = iyj2eB(n + l)I n+ltk (x). 
Combining eqs. (§) and @ yields 

= I n ,k( x )> E 2 = m 2 + 2eBn. 

The function g K (x) satisfies 



9k{x) 



E + m 



fn(x), 



implying that 



9k(x) = -iV2eBnI n _ hk (x). 
The normalized eigenfunctions become 



WaCM) = exp^Ent + iky)^—^ — Ttv ^ 



2E„ 



I n -i,k{x) 



(10) 

(11) 
(12) 

(13) 



where n = 0,1,2,..., E n = y/m 2 + 2eBn and ip^(x.,t) are positive and negative 
energy solutions, respectively. Note that ipQ~ k (x, t) = and that we have defined 
I-i,k( x ) = 0. The spectrum is therefore asymmetric and this asymmetry is inti- 
mately related to the induced vacuum charge, as will be shown in section 4. In fig. [1] 
a) we have shown the spectrum for m > and in fig. [I] b) for m < 0. 



The field may now be expanded in the complete set of eigenmodes 

CO - 

¥(x,t) = £ / dk \b n ^k\*,t)+d: jk i; { - k \x,t) 



n=0 



(14) 



Quantization is carried out in the usual way by promoting the Fourier coefficients 
to operators satisfying 

{b n ,k,b n > k ,} = S n)n i8k }k /, {d n ,k,dl l , k ,} = 5n tn '5k t k', (15) 
and all other anti-commutators being zero. 
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3 Fermion Propagators and the Effective Action 



In the previous section we solved the Dirac equation and with the wave functions 
at hand, we can construct the propagator. From the trace of the propagator the 
effective action to one-loop order is calculated. 



The fermion propagator. The fermion propagator in vacuum is 

iS F (x',x) = (0| 7 1 [tf(x\t / )*(x,t)] |0), 
where T denotes time ordering. By use of the expansion fll4|) one finds 



(16) 



iSf^x' , x) 



~ rdk 

hJ ^ 



(17) 

After some purely algebraic manipulations and using the integral representation of 
the step function, we obtain 



S F (x', x) ab = -— ^ V / dkduj En ^ m exp [ -iu{t' -t) + ik(y' - y) ] 

^ 2 ^Q J 2E n 



uj 2 -EI + xe 
(18) 



S ab (n,u, k). 



Here S a b(n : uj, k) is the matrix 

/ I nik (x ')I n)k (x) y/ 2eBnI n)k (x 1 ) I n _ 1>k (x) 

V E ^f^V2eBnI n _ hk (x')I ntk (x) I n -i,k{x')In-i,k{ x ) 



(19) 



At finite temperature and chemical potential we write the thermal propagator as 
(see ref. [[Hj for details) 

(S F (x',x))^ = S F (x',x) + S^(x',x). (20) 



The thermal part of the propagator is 



F V A 5 



X , X) 



£ 

n=0 



dk 
2^ 



(21) 



where 



(+)/ 



1 



exp + //) + ! 



(22) 



exp (5{u — fx) + 1 

This may be rewritten as 
°° f dk 

)=iV / — expzfc(y'-y)expiu;(t / -t)/ F (c ( ;)5(u; 2 -^-2£:)5'f f) (n,u;,fc). 

(23) 
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Here S^ b (n,uj, k) is the matrix 

/ (lj + m)I ntk (x')I n)k (x) 



In-l,k(x')InA X ) 



In,k( x ')In-l,k( x ) ( U ~ rn)In-l,k(x')In-l,k(%) 



(24) 



As noted in ref. [TIJ, one is not restricted to use equilibrium distributions in this 
approach. Single particle non-equilibrium distributions may be more appropriate if 
e.g. an electric field has driven the system out of equlibrium. 

The effective action. The generating functional for fermionic Green's functions in 
an external magnetic field may be written as a path integral: 



Dip T>ip exp 



d 3 x {--F iu/ F ttv + m)ip -fjip + tprj) 



(25) 

The functional integral describes the interaction of fermions with a classical elec- 
tromagnetic field. It includes the effects of all virtual electron-positron pairs, but 
virtual photons are not present, which means that we are considering the weak cou- 
pling limit. 

The fermion field can be integrated over since the functional integral is Gaussian: 

1 



Z(r], t), A^) = det[i(ip— m) ] exp 



-F, U F^+ I d"yfj(x)S F (x,y) V (y) 



(26) 

Taking the logarithm of Z(rj, rj, A^) with vanishing sources gives the effective action 



S, 



eff 







f d 3 x 





iTr log [i(ip— m) 



(27) 



where we have written logdet = Tr log by the use of a complete orthogonal basis. 

m yields 

itrSp(x, x). 



Differentiating eq. (27) with respect to m yields 



dm 



(2? 



The trace is now over spinor indices only. By calculating the trace of the propagator 
and integrating this expression with respect to m thus yields the one-loop contribu- 
tion to the effective action. This method has been previously applied by Elmfors et 
al ]l(J in 3+1 dimensions. The above equation may readily be generalized to finite 
temperature, where we separate the vacuum contribution in the effective action 



C = Co + Cx + C^ 1 = C + C 



eff 



where Cr, is the tree level contribution, and 



dC 



eff 



dm 



i tr 



Sf(x,x) + Sf(x,x) 
6 



(29) 



(30) 



Using eqs. (|18|) and (|19|) a straightforward calculation gives for the vacuum contri- 
bution 



trS F (x,x) = jr J ^ dk ^ + ie [m (ll >k (x) + i^_ 1[fc (x)) + u (i^ |fc (x) - i, 



1 y g?/cgL> 2 



n—l 



ieB 



Integrating this expression with respect to m yields 

C x = -— y Vm* + 2eBn - (32) 
2tt n=1 4vr 

The divergence may be sidestepped by using the integral representation of the 
gamma function [13| and subtract a constant to make Ci vanish for B = 0, 

1 r°° ds 

C\ = — 3- / — exp(— m 2 s) [eBs coth(eBs) — 1]. (33) 

8-7T2 JO S 2 

This result calls for a few comments. Firstly, one notes that the induced Chern- 
Simons term in the effective action vanishes. This is not in conflict with the result 
of Redlich; rather it is caused by our choice of gauge, A = 0. Insertion of A = into 
the general expression for the Chern-Simons term obtained by Redlich, will make 
it vanish. Thus, our results agree. Secondly, one can show that the Chern-Simons 
term is gauge-invariant up to surface terms and using another gauge the Chern- 
Simons term needs not vanish. The choice Aq ^ gives rise to a Chern-Simons 
term AoBme 2 / (\m\8ir 2 ) in this approach and the result is again in accordance with 
the general expression. It is always satisfactory to see that identical results can be 
obtained by entirely different methods. 

The finite temperature part of the effective action is calculated analogously using 
the thermal part of the propagator (|2~ID, 



rp r> oo 

= — E [ M 1 + ex P ~P(E n -//)] + log[ 1 + exp -f3{E n + //)]] (34) 

/7T L J 



2tt n=1 



TeB 

— — log [ 1 + exp -f3(m -/d)]. 

Z7l 



Letting B — > it can be shown that one obtains the free energy of a gas of non- 
interacting electrons and positrons: 



= — I EdE 

Z7T JO 

where E = \/m 2 + k 2 . 



T 



log[ 1 + exp —f3(E -fi)] + log[ 1 + exp -(3(E + n) 



(35) 



In the following we restrict ourselves to the case \x > 0. Analogous results can 
be obtained for \i < 0. 



In the zero temperature limit of C^ ,fl one gets 

271 n=0 



(36) 



where the prime indicates that the sum is restricted to integers less than (/i 2 
m 2 )/2eB. Similarly, one may derive the charge density at T = 0: 



P 



dCP* eBr fi 2 -m 



djj, 



2tt 



Int(- 



2eB 



+ 1 



ji> m, 



(37) 



in accordance with the result of Zeitlin \TM . Notice that at T = the charge density 
is equal to the particle density, since no antiparticles are present, as can be seen by 
inspection of eqs. (|36l) and (|37|). 



4 Magnetization and the de Haas- van Alphen Ef- 
fect 



In this section we study the physical content of the effective action which was ob- 
tained in the previous section. In particular we investigate a few limits to check the 
consistency of our calculations. 



The magnetization is defined by 



M 



dC 



eff 



dB 



The vacuum contribution to the magnetization is obtained from eq. 



Ml 



87T2 Jo 



°° exp {—m s , 
as— 



5 
S'2 



es coth(e-Bs) — 



e 2 B.r 



sinh 2 (e73s) 



(3* 



(39) 
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For the thermal part of the magnetization we find 



M 13 " 



Te 



— Y] \ log[ 1 + exp -(3{E n -n)\+ log[ 1 + exp -(3(E n + /i) 

/7T L 



2tt . 

n=l 

Te 

H log [ 1 + exp —(3 (ra — jti) ] 



2tt 

e 2 5 ~ 



n 



1 



exp (3(E n — li) + 1 exp /3(i? n + /i) + 1 



(40) 



Magnetization at zero temperature. In the zero temperature limit eq. (|40|) reduces 
to 



2tt 



E 

n=0 



eBn 
E„ 



(41) 



where the sum again is restricted to integers less than (/x 2 — m 2 )/2eB. 

m 2 m 2 ) the vacuum contribution 
e 2 B 



In the weak 5-field limit [eB fi 2 
becomes 

e 2 B f°° , exp(-m 2 s) 



127T 771 



(42) 



This agrees with the results of refs. [14,15]. The contribution from real thermal 
particles is obtained by rewriting the square root in C 13 '^, using the integral repre- 
sentation of the gamma function and treating Int([i 2 — m 2 )/2eB as a continuous 
variable in the limit B — > 0. The result is 



M 13 " 



4tt 



(fx 



771 



(43) 



Some comments are in order. It is perhaps somewhat surprising that the magnetiza- 
tion is non-zero in this limit. One should, however, bear in mind that the sign of m 
uniquely determines the spin of the particles (and antiparticles), implying that the 
system under investigation consists entirely of either spin up or spin down particles. 
This is not the case in 3+ld, where the representations characterized by the sign of 
771 are equivalent. It would therefore be natural to consider a system consisting of an 
equal number of spin up and spin down particles, which then amounts to sum over 
±m []. By doing so one finds a vanishing magnetization as B goes to zero, exactly 
as in 3+1 dimensions. Note that this result, of course, can be obtained by using 
four component Dirac spinors, since the four dimensional representation reduces to 
the two inequivalent two dimensional ones. 



1 Onc should use the same chemical potential irrespective of the sign of m in order to ensure 
the same charge density for spin up and spin down particles. 
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In order to get the strong field limit ( eB 3> /i 2 , m 2 ) of the vacuum contribution, we 
scale out eB and take eB — > oo in the remainder. This gives 

d oc (eB)t => M x oc eiv^B. (44) 

Vacuum effects contribute to the magnetization proportional to the square root of 
B. This should be compared with the corresponding result in 3+ld, where the 
magnetization goes like Blog(^) flQfl . For real thermal particles only the lowest 
Landau level contributes in the strong field limit, as can be seen directly from eq. 
([4l"l) 0. The thermal part of the magnetization is then M^ ,fl = e(/i — m)/27r. For 
large S-fields all particles are in the lowest Landau level and [i = m. Hence, the 
contribution to the magnetization from real thermal particles vanishes. This can be 
understood from the following physical argument: The energy of the single particle 
ground state is independent of the external field (Eq = m), so increasing B cannot 
lead to an increase in £^ ,At , when the charge density (and therefore the particle 
density) is held constant. 

High Temperature Limit (T 2 ^> m 2 ^> eB,fi = 0). The high temperature limit 
is rather trivial. From a physical point of view, one expects that Cr ,tl approaches 
the thermodynamic potential of a gas of non-interacting particles of mass m. Indeed 
we obtain 

mT 2 T 3 
C 0,» = -_L* 2 [ l + exp(-/3m) ] + — Li 3 { 1 + exp(-/?m) ], (45) 
2ir 2rc 



which, of course, coincides with eq. (|35|) when /i = 0. Here Li n (l + x) is the 
polylogarithmic function of order n: 

oo k 

U n (l + x) = Y,-- (46) 
fc=i K 

Eq. (|4"5| ) may be obtained by carefully replacing the sum by an integral. 

Magnetization at finite temperature. In fig. ||] we have displayed the total mag- 
netization as a function of the external magnetic field for different values of the 
temperature (/x/m = 1/1.5, T/m =1/150 solid line, 1/50 dashed line, 1/5 dotted 
line). Fig. ^ is a magnification of fig. ^| in the oscillatory region. The de Haas- van 
Alphen oscillations are seen to be present for low temperatures and weak magnetic 
fields. Finally, we have summed the magnetizations for ±m. The resulting curve 
(T/m= 1/100) is shown in fig. |] and reveals that the magnetization goes to zero in 
the limit B -> 0. 



2 This is also the case in 3+ld where the energy of the lowest Landau level is independent of B. 
Eq = \J to 2 + kl , where k z is the z-componcnt of the momentum. 
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Induced vacuum charges and currents. Finally, we calculate the vacuum expec- 
tation value of the induced charge and current densities. Such calculations have 
been carried out in other contexts, e.g. in connection with magnetic flux strings 
(see ref. ||). We shall employ the most commonly used definition of the current 
operator which can be shown to measure the spectral asymmetry relative to the 
spectrum of free Dirac particles. 



f (x) 



(47) 



The induced charge density is given by 



(f(x)) = (p(x)) 



--E 



dksign(E n ) \ ^ k (x) 



(4f 



Using the complete set of eigenmodes as given by eq. (|i"3|) , a straightforward calcu- 
lation yields 

, . . . m e 2 B . . 

(p(x)) = ri — . (49) 
| m\ 47r 

Eq. fl4*9|) is simply the Chern-Simons relation and our calculations are thus in com- 
plete agreement with the result of Fumita and Shizuya. This result has the following 
physical interpretation: As we turn the magnetic field on, an unpaired energy level 
E = m emerges (in the case m > 0). The number of spin down negative energy 
electrons therefore gets reduced relative to the free case by a factor |-, which is 
the degeneracy per unit area. This can be interpreted as the appearance of spin up 
positrons and results in a positive charge density. For m < a similar argument 
applies. 

A similar calculation for ( j (x) ) reveals that the induced current vanishes. This 
result should come as no surprise due to translational symmetry of the system. A 
non-vanishing vacuum current would arise in the presence of an external electric 
field and is then attributed to the drift of the induced vacuum charge. 



5 Conclusions 



In this paper we have calculated the effective action for fermions moving in a plane 
with a constant magnetic field orthogonal to the plane. We have derived the mag- 
netization from the effective action and have shown that the system exhibits the 
well-known de Haas-van Alphen oscillations at small temperatures and low values 
of the magnetic field. 
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Finally, it would be of some interest to extend the present work. Firstly, one 
should examine improvements to our results by considering the corrections to the 
self-energy of the electrons in the presence of fields and a thermal heat bath. Sec- 
ondly, one could treat the highly non-trivial problem of fermions in slowly varying 
electric and magnetic fields. 

The authors would like to thank M. Burgess for useful comments and sugges- 
tions. 
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FIGURE CAPTIONS: 

Figure 1: The energy spectra of Dirac fermions in the presence of a magnetic field, 
a) m > and b) m < 0. 

Figure 2: The magnetization in units of em as a function of B in units of m 2 /e for 
different values of temperature. \x = 1.5m. 

Figure 3: A magnification of the oscillatory region in fig. 2. 

Figure 4: The magnetization for a system consisting of an equal number of spin up 
and spin down particles. 
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